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The electric current conservation in a two-dimensional quantum wire under a 
time dependent field is investigated. Such a conservation is obtained as the global 
density of states contribution to the emittance is balanced by the contribution 
due to the internal charge response inside the sample. However when the global 
partial density of states is approximately calculated using scattering matrix only, 
correction terms are needed to obtain precise current conservation. We have 
derived these corrections analytically using a specific two-dimensional system. 
We found that when the incident energy E is near the first subband, our result 
reduces to the one-dimensional result. As E approaches to the n-th subband with 
n > 1, the correction term diverges. This explains the systematic deviation to 
precise current conservation observed in a previous numerical calculation. 
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1 Introduction 



The dynamic conductance of a quantum coherent mesoscopic system under a time 
dependent external field is the subject of recent interests [[I], 0, |3|, f|, In con- 
trast to dc-transport where the internal potential distribution inside the sample 
does not appear explicitly, the AC-response depends sensitively on the internal 
potential distribution. This internal potential is due to the charge distribution 
generated by the applied AC-field at the leads and it has to be determined self- 
consistently|l]]. So far there are two approaches to the coherent AC-transport 
problem. One is to derive a formal linear response to a given potential distribu- 
tion in the sample @. The difficulty with such an approach is that the poten- 
tial distribution is not known a priori. Another approach is to investigate the 
AC-response to an external perturbation which prescribes the potentials in the 
reservoirs only]?], |J. The external potentials effectively determine the chemical 
potential of the reservoirs and the potential distribution in the conductor must 
be considered a part of the response which is to be calculated self-consistently. 
In this approach, Biittiker and his coworkers |], |§ have formulated a current con- 
serving formalism for the low frequency admittance of mesoscopic conductors. 

In the theory of Biittiker, Pretre and Thomas|l|, it is necessary to consider the 
Coulomb interactions between the many charges inside the sample, in order to 
preserve the current conservation. For a multi-probe conductor the low frequency 
admittance is found to have the form|| |9| G a p(u) = G a p(0) — iujE a/3 + 0(u; 2 ), 
where G a /3(0) is the dc- conductance, E a p is the emittance||, and a (or (3) labels 
the probe. The emittance E a/3 describes the current response at probe a due 
to a variation of the electro-chemical potential at probe (3 to leading order with 
respect to frequency u. It can be written as|| || E a p = dN a p/dE — D a p, where 
the term dN a/3 /dE is the global partial density of states (GPDOS) flTOfl which is 
related to the scattering matrix. It describes the density of states of carriers 
injected in probe j3 reaching probe a and is a purely kinetic term. The term D a p 
is due to the Coulomb interaction of electrons inside the sample and is a term 
of capacitive nature. D a $ can be computed from the local density of states |], Hi 
which is related to the electron dwell times. Electric current conservation, namely 
J2a G a p{uj) = 0, means that J2 a E a p = or equivalent ly|], [TTH 

dN p _^dNap n T dj/3 m 

~dE = \^E-\ Da!i -~ (1) 

where dNp/dE is the DOS and is the dwell time for particles coming from 
the probe (3. Clearly the current conservation is established since one realizes 



2 



that J2 Q dN Q/3 /dE is the physical quantity called injectance which is identical^]] 

to E a D a/3 . 

Applying the above formalism to mesoscopic conductors, one needs to com- 
pute various physical quantities such as the partial density of states. These 
quantities have vivid physical meaning|| but are not easy to obtain exactly. 
For a large system, the GPDOS can be expressed approximately in terms of the 
energy derivative of the scattering matrix elements ||12||: 



dN, 



a/3 



dE 




(2) 



Because for a given system one may be able to obtain the scattering matrix, Eq. (0) 
thus provides a practical means of computing the GPDOS. On the other hand, 
in order to obtain current conservation precisely, a correction should be added 
to Eq. (D) which can be neglected for large systems and large energies ||T0|, |i3| . 
For one-dimensional systems, such a correction has been derived by Gasparian 



et. al. [13fl which contains the reflection amplitude divided by the energy, 



dNa _ dNa j (3) 

dE dE l 4nE J ' V ; 

where dN a /dE = J^pdN^/dE which is computed from Eq. (H). 

We have recently applied the above current conserving formalism to a two- 
dimensional mesoscopic conductor in the shape of a T-junction^. To the best of 
our knowledge, it was the first 2D calculation with first principles. Among other 
things, an interesting and we believe useful discovery was that Eq.(U) turned out 
to be inaccurate in 2D. First of all, energy E in the second term on the right 
hand side of Eq.@ has to be replaced by the longitudinal part of the incident 
energy. Even with this change, there were small but systematic deviations to 
precise current conservation when the energy is approaching the second subband. 
In fact it was found that the DOS dN a /dE as defined above diverges near the 
onset of the second subband and this led to the observed systematic deviations^. 

We are not aware of any 2D theory to account for the correction term which 
should appear in Eq. (|3]). The purpose of this paper is to investigate such 
correction terms in two dimensions. This not only provides further theoretical 
insights to the problem of AC-transport, but is also very helpful from a practical 
application point of view. From our own experience, numerical AC-transport 
calculations can be quite tricky and being able to obtain precise electric current 
conservation often serves as a very stringent check to numerical results. To this 
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purpose, we have considered the simplest two-dimensional model which is a 8- 
potential inside a quasi- ID ballistic conductor Jl^]. Since quantum scattering in 
this system leads to mode mixing which is the basic feature of a two-dimensional 
system, it provides answers to our 2D problem. The advantage of this system is 
that it can be solved exactly. We have thus derived analytically the correction 
term. In particular we found that when the incident energy E is within the first 
subband, our result essentially reduces to the one-dimensional result Eq. (|3]). As 
E is increased to approach the n-th subband edge with n > 1, the correction 
term diverges. This explains the systematic deviation observed in our previous 
numerical calculation p. 

The paper is organized as the following. In the next section we present the 
solution of the 2D scattering problem and derive the correction term. Section III 
contains our numerical tests of the analytical formula. The last section serves as 
a summary. 



2 Model and results 

Figure 1 shows the system where a ^-potential is confined inside a quasi- ID wire 
with width a. We assume, for simplicity of the calculation, that the boundaries 
of the ballistic conductor are hard walls, i.e. the potential V = oo. Inside the 
conductor, the potential is zero except that a 8 function potential V(x, y) = 
■y8(x)8(y — yo) is placed at r = (0,yo). The scattering region x\ < x < X2 is 
assumed to be symmetric with x 2 = = L/2. From now on we set h — 1 and 
m = 1/2 to fix our units. 

To compute the transmission and reflection amplitudes thus the scattering 
matrix, a mode matching method fi~5|l is employed. The electron wave functions 
are written as follows. For region I (see figure 1): 

* / = E Xn(y) (a n e iknX + b n e~ iknX ) , 

n 

where Xn{y) is the wave function of the n-th subband along y-direction; a n is the 
incoming wave amplitude and taken as an input parameter; h n is the reflection 
amplitude; and = E — (nn/a) 2 is the longitudinal momentum for the n-th 
mode. Note that for electron traveling in the first subband, k n with n > 1 is 
purely imaginary. Similarly for region II: 

= E Xn(y) (c n e iknX + d n e~ ik " x ) , 
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where c n is transmission amplitude and d n is set to zero in our calculation. After 
matching the boundary conditions at x — 0, we obtain 

O'n b n c n 

and 

ik n c n ik n [a n b n ) ^ ] I" ^ Jim (o m + & m ) , 
where T nm = 7*£G/o)Xm(yo)- Eliminating c n , we have 

e = P6 , (4) 

— * 

where e n = -Y, m ^nma m and P nm = T„ m - 2ik n 5 nm . To find b we need to 
compute P^ 1 . Introducing a new matrix P = I + M with M„ m = ir nm / (2k m ) so 
that P nm (—2ik m ) = P nm . Expanding P -1 in powers of M, we have 

P" 1 = — - — I ■ M ■ M 2 - M 3 ... 
I + M 

Since r nm T m i = T n iT mm , we find that M 2 = (a-l)M where a = 1+z En E„ n / (2k n ), 
from which we have P _1 = 1 — M/ct. Finally, we obtain the matrix elements 

{P )nm — „, ($nm ~ ) ■ (5) 

We shall specialize to consider the incident electron being in the first subband: 
a n = <W Using Eq.fl^D the reflection and transmission amplitudes are 

b n = ^(P'^nmem = ^ ^ , (6) 

c n = 5 n i + b n . (7) 

For our system the scattering matrix elements s a p are given by Sn = 61 exp(ikiL) 
and S12 = ci exp(ikiL). The approximate DOS becomes, using Eq. (0), 

dA^ 1 ^ ( + cfe Q/3 ds a/3 



dE Aui y \ ap dE dE 

~ Im \ J -iEjMr • ( 8 ) 



47rA;i \47ik 2 J Am ^ %k\k n 

To derive this expression we have used a relation 26^ + 1 = a /a* which follows 
directly from the unitary condition of the scattering matrix. Next we compute 
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the dwell time and hence the precise DOS (as opposed to the approximate DOS 
ofEq. (§): 



T di = / \^i\ 2 dxdy + / \^fjj\ 2 dxdy 
Ji J ii 



T p ik\L _ i p ik n L _ i 

From Eqs. (HD , (§) and (H), we arrive at the following central result of this work, 
dNa _ dN» T r Saa i , J_ V 1 &W I 2 /-, n x 

Hence we found that for this 2D system, there are two correction terms to the 
DOS. Clearly the first correction term, i.e. the 2nd term on the right hand side 
of Eq. (|T0|) , is generic, as it can be written in terms of the scattering matrix 
element. This term is similar to the corresponding term in Eq. @ of the ID 
case, except that the total energy E in Eq. @ is now replaced by the transport 
energy k\. In fact this them has been guessed in our earlier work||. There is a 
second correction term (the 3rd term of Eq. (|T0|) ) which comes solely due to mode 



mixing in our 2D system, and understandably it does not exist in ID cases [ 13]. 

For small incident energies, i.e. as k± goes to zero, |6 n | 2 — > kf for n > 
1. Therefore the second correction term of ( |10D is actually negligible at small 
energies. Indeed, this is the case in our earlier numerical calculations |J where 
current conservation was very well satisfied at low energies using Eq.(|3|). However, 
as energy is approaching the n-th subband edge, for small k n — > with n > 1, 
|6 n | 2 remains finite. Hence according to Eq. ( |Tol) the second correction term 
diverges at these higher subband edges. This explains the observation of our 
calculation || where systematic numerical errors exist in current conservation near 
the 2nd subband edge. For energies within the first subband, as mentioned above 
k n are all pure imaginary numbers with n > 1. Hence with large system size L, 
the factor exp(ik n L) is very small as long as k n ^ 0. However we emphasis that 
the second correction term becomes dominant very near each subband edge thus 
must be included in order to obtain precise current conservation. 

Finally we note that Eq.fllOD is not coordinate independent, so care must be 
taken when using it. For instance, if we choose X\ as the origin in figure 1, the 
factor exp(ik n L) in the last term of Eq.(|H]) will be canceled due to the coordinate 
shift while the second term of Eq.flTUD remains the same. In this sense, the new 
correction term is not generic and must be computed case by case for 2D systems. 
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3 Numerical test 



To gain further intuitive impression of the AC-transport, and in particular to 



check our analytical formula Eq.(lO), we shall first present direct numerical cal- 
culations of the admittance for the quantum wire system studied in the last 
section (Fig.(l)). Obviously since this problem was solved exactly above, agree- 
ment is obtained with Eq. flTTj) . We shall then study the validity of Eq.flTOD using 
another more complicated 2D conductor in the shape of a T-junction (see be- 
low). Indeed, although Eq.fllOD was derived using a specific example of Fig.(l), 
it dramatically improves the current conservation near the second subband edge 
for the T-junction as well. 

In order to compute the admittance, we have to know D a ^ which is related to 
the dwell timefl], For a metallic conductor, it is appropriate to use the Thomas- 
Fermi approximation. Under such an approximation D a ^ is given by|l], || 



,, (dn(a,r)/dE)(dn(r,0)/dE) 



where the local density of states dn(r, [3)/dE is the injectivity which measures 
the additional local charge density brought into the sample at point r by the 
oscillating chemical potential at probe (3. The injectivity can be expressed as[[TJ 

dnjrj) _ |^ n (r)| 2 

dE ^ 27rv Pn ' 1 ] 

where vp n is the velocity of carriers at the Fermi energy at mode n in probe (3. 
dn(a, r)/dE is called the emissivity which describes the local density of states of 
carriers at point r which are emitted by the conductor at probe a. It is defined 

as 

dE ~ Amy T Saf3 eSU(r) e8U{r) SaP 

It has been shown |J that in the absence of magnetic field the injectivity is equal 
to the emissivity. Using Eqs.(||), (|TTD and (|12|), we can calculate the emittance. 

Specifically, for the system of Fig.(l) we consider incident electron coming 
from probe 1 and set a = L = 1, yo — 0.3, and 7 = — 1. In Fig. (2), we plot 
the global DOS together with the transmission coefficient T. As expected, the 
transmission coefficient T(E) (solid line) has large values for almost all energies 
E except at a special energy E r where we have complete reflection (reflection 
coefficient R(E r ) = 1) due to the resonant state. This can also be seen from the 



(13) 
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behavior of the global partial DOS for reflection dNu/dE (dotted line) which 
peaks when T(E = E r ) = 0. On the other hand, dN 2 \/dE (dashed line), which is 
the global partial DOS for transmission, takes minimum value at E = E r . This 
behavior is consistent with that of a ID system made of a symmetric scatterer ||10|| 
where one has dN n /dE ~ RdN/dE and dN 2l /dE ~ TdN/dE. In Fig. (3), the 
quantities Dn (solid line) and Dyi (dotted line) are shown. Both curves reach 
maximum values near the resonant point E r , which is expected since D a p are pro- 
portional to the dwell time or the DOS. The emittance E a p is plotted in Fig. (4). 
Both E\\ (solid line) and E\ 2 (dotted line) reach extremal values at the resonant 
point. The system responds differently for different energy, either capacitively 
when E\\ = —E\ 2 > 0, or inductively otherwise. From Fig. (4), we observe that 
near the resonance E u and E 12 respond capacitively while E\ 2 is inductive away 
from this resonance energy. This behavior, namely being capacitive when at the 
T « resonance, is the same as that observed in the 2D T-junction^. On the 
other hand for an ID tunneling system[|l|] the response is inductive at its reso- 
nance. But in that case the resonance is marked by transmission coefficient being 
near unity. Finally, to confirm electric current conservation, essentially the two 
curves of Fig. (4) must add to zero. Clearly these curves do not cancel each other 
as the figure shows, exactly due to the approximate nature of the partial density 
of states as obtained using Eq. (0). After including the two corrections to DOS 
as derived in Eq. (|10D , however, we did obtain perfect current conservation for 
the whole energy range. This is not surprising since after all (|H]) is an exact 
result for this quantum system. 

Our main result Eq.(|T0|) is derived using a specific simple example shown in 
Fig.(l). There seems no special reason for Eq.(|K]) to apply to other 2D systems, 
since the form of the new correction term is given by the amplitudes of the non- 
propagating modes inside the scattering junction (as oppose to the more general 
scattering matrix elements), and these evanescent amplitudes probably depend on 
the scatterer in some fashion. In this sense it is unfortunate that a more general 
form was not obtained. However since the new correction term does explain, 
qualitatively, the observed discrepancy of using Eq. @ as discussed above, it 
is tempting to test it using the more complicated 2D system of the T-junction 
studied previously 0. As the T-junction has been reviewed and studied by many 
authors |16|, [5] at various contexts, here we shall not present the details for its 
calculation. For this purpose, we have checked the current conservation of the 
T-junction using Eq.(^|) and compared with the result obtained using Eq. (|T0l) . In 
Fig. (5), we have plotted the DOS dN x jdE given by Eq.(|) (dotted line) and bypTf 
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Eq. (|10D (solid line), and the dwell time r^i/ln (dashed line). Although the result 
Eq. fllCf ) is model dependent, we observe that the agreement is clearly better. This 
suggests that the new correction term does capture the essential ingredient of the 
correction, although it is not completely universal as the evanescent amplitudes 
depend on the peculiarities of a 2D system in some weak way, leading to the small 
remaining difference. 

4 Summary 

In summary, we have investigated the electric current conservation in a two- 
dimensional ballistic conductor under a time dependent field. Similar to that of 
the ID case, we found that in order to obtain precise current conservation, certain 
corrections to the density of states as obtained approximately from the scattering 
matrix must be included. We have derived these corrections analytically for a 
specific two-dimensional system and found that there are two correction terms. 
One of the correction term has the same form as that of the ID case, while the 
second correction term is purely due to mode mixing characteristic of 2D quantum 
scattering. In particular, when the incident energy E is within the first subband, 
our result essentially reduces to the one-dimensional result if E is not too high. On 
the other hand as E approaches to the n-th subband with n > 1, the correction 
term diverges at the subband edges. Hence in 2D the mode mixing leads to 
important changes in the global density of states and must be included if precise 
electric current conservation is desired. Finally, the new correction term found 
here provides a qualitative explanation for the small but systematic deviation to 
precise current conservation observed in our previous numerical calculations [[|] on 
a 2D quantum wire in the shape of the T-junction. Indeed, our numerical test 
has produced better agreement when the new formula derived here is used. 
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Figure Captions 

Figure 1. Schematic plot of the quantum wire system: a 5 potential 75(r* — r~o) is 
confined inside a quasi-lD quantum wire, with fo = (0, yo). The wire width 
is a. The scattering region is between x\ and where x 2 = — #1 = L/2. 
In our calculations, the parameters are set to L = a = 1, yo = 0.3, and 
7 = -1.0. 

Figure 2. The global partial density of states and the transmission coefficient as func- 
tions of electron energy E. Solid line: transmission coefficient T; dotted 
line: dNu/dE; dashed line: dN2\/dE. Unit of energy is h 2 /2ma 2 . 

Figure 3. The current response to the internal potential, D a/3 , as a function of energy 
E. Solid line: dotted line: D 2 i- 

Figure 4. The dynamic part of the admittance, E a p = dN a p/dE — D a p as a function 
of energy. 

Figure 5. A numerical check of the electric current conservation, Eq. ([1]), for the T- 
junction studied in Ref. |J. Solid line: dNi/dE as obtained by Eq. (|10|); 
dotted line: r^/h = J2 a D a i/h. Agreement of the two curves indicate the 
conservation. The remaining small differences at high end of the energy 
between the two curves indicates that the new correction term in Eq. ( |T0"D 
has a weak non-universal dependence on the 2D system shapes. 
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